Scintillator Non-Proportionality and Gamma Quenching in CaWC>4 
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7 Deceased 

We measure and explain scintillator non-proportionality and gamma quenching of CaWC>4 at low 
energies and low temperatures. Phonons that are created following an interaction in the scintillating 
crystal at temperatures of ~ 15 mK are used for a calorimetric measurement of the deposited 
energy, and the scintillation light is measured with a separate cryogenic light detector. Making use 
of radioactivity intrinsic to the scintillating crystal, the scintillator non-proportionality is mapped 
out to electron energies < 5keV. The observed behavior is in agreement with a simple model 
based on Birks' law and the stopping power dE/dx for electrons. We find for Birks' constant 
ks = (18.5 ± 0.7) nm/keV in CaWCU. Gamma lines allow a measurement of the reduced light yield 
of photons with respect to electrons, as expected in the presence of scintillator non-proportionality. 
In particular, we show that gamma-induced events in CaWCU give only about 90% of the light yield 
of electrons, at energies between 40keV and 80keV. 

PACS numbers: 29.40.Mc,29.40.Vj 
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I. INTRODUCTION 

Scintillation is among the oldest and most classic of 
particle detector technologies. For photons and electrons 
with energies below ~ 100 keV, it appears that the light 
output from the scintillator is not simply proportional to 
the absorbed energy [1]. Recently, this so-called scintil- 
lator non-proportionality is receiving renewed attention 
thanks to the utilization of scintillators at these low ener- 
gies in the search for dark matter [2] . For ions and other 
particles producing heavily ionizing tracks on the one 
hand, this effect can be understood as a saturation effect 
that occurs for the light produced per unit length dL/dx 
when an increase in local energy deposition dE/dx does 
not lead to more light. This phenomenon is encoded in 
Birks' law [3] 

dL _ AdE/dx 

dx ~ l + k B dE/dx' ^ ' 

A and ks are empirical constants that depend on the 
material and have to be determined from the data. At 
high energy loss dE/dx, the light output is no longer 
proportional to the energy loss, due to the term kgdE/dx 
becoming larger than unity. Eventually, dL/dx becomes 
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independent of dE/dx and the total light from an event 
depends only on the track length [4] . 

Electrons on the other hand are lightly ionizing, but for 
low-energy electrons, there is in fact an increase in local 
energy loss. As an example, figure 1 shows dE/dx (E) 
for electrons in CaW04, the material we shall study here. 
One might hence anticipate some effect from the non- 
proportionality in equation 1 also for low-energy elec- 
trons, where dE/dx for energies near the end of the track 
(i.e. below ~ lOkeV) is significantly increased. 
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FIG. 1: Energy dependence of d_E/d:r for electrons in CaW04, 
from the estar database [5]. 

Evidently, the most direct way to measure the response 
of a scintillator to electrons is to study its light output 
from single monoenergetic electrons, over a range of low 
energies. Since there are no natural electron sources of 
this type, the Compton-coincidence technique has been 
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developed [6]. In this technique, one uses a monoener- 
getic gamma with relatively high energy, which under- 
goes a Compton scattering in the material under study. 
The scattered photon is observed in a separate detec- 
tor. Knowledge of the scattered photon's angle leads, via 
Compton scattering kinematics, to a recoil electron of 
known energy in the material under study. This method 
has been employed successfully to map out the response 
of various scintillators to low-energy electrons (see [7] and 
references therein). 

Another method to measure the electron response of a 
scintillator uses the monochromatic X-ray radiation from 
a synchrotron source. The scintillation light yield of an 
excitation above the characteristic K-edgc of the clement 
is measured. When the energy independent response 
from the cascade event that is due to the relaxation of 
holes in the K-shell is subtracted, it is possible to derive 
the response due to the K-shell photoelectron alone [8, 9]. 

In this note, we would like to draw attention to the 
possibility of a new method, which has arisen through 
the investigation of cryogenic scintillating detectors in 
the search for dark matter. A two-channel readout al- 
lows both the energy and the light yield of an event to 
be determined. A large sample of low energy electron 
events with both known energy and light yield, acquired 
as a background in the search for dark matter, allows a 
precise measurement of the scintillator response and its 
non-linearity. 

The fundamental processes of excitation of elec- 
trons in deep shells of atoms are temperature inde- 
pendent. Indeed, it has already been shown that the 
non-proportionality of self-activated scintillators such as 
CaWC>4 or BGO exhibit no temperature dependence [10, 
11] even thought their light yield changes significantly 
with temperature [12, 13]. Hence, measuring at mil- 
likelvin temperatures allows one to infer the scintillator 
non-proportionality also at room temperature. 



II. SETUP 

The CRESST-II experiment [14] is built to search for 
and detect a new particle species, so-called weakly inter- 
acting massive particles (WIMPs), which could consti- 
tute the dark matter known to exist in the universe [15]. 
The expected interaction rates are less than one event 
per kilogram of target material, per year of observation 
and per keV in the energy range between ~ lOkeV and 
~ 40keV [14]. This requires a substantial reduction of 
backgrounds from ambient radioactivity. To this end, 
the CRESST-II experiment is located in the Laboratori 
Nazionali del Gran Sasso under an average rock overbur- 
den of 1400 m to shield the experiment from cosmic rays, 
and additionally, the cryostat with the target materials 
is surrounded by thicknesses of 45 cm of polyethylene, 
20 cm of lead, and 14cm of copper. 

The remaining activity is mostly due to beta emitters 
intrinsic to the employed target materials [16]. To distin- 



guish the expected WIMP-induced nuclear-recoil signal 
from this remaining background, an active discrimination 
technique is used. In CRESST-II, we use scintillating 
tungstates (CaW04, ZnWO^) as target materials. These 
crystals are shaped as cylinders with a diameter of 4 cm 
and similar height, thus weighing about 300 g (CaWO.4) 
or 400 g (ZnWO^) each. They are cooled to temperatures 
of ~ 15 mK and operated as cryogenic calorimeters [17]: 
non-thermal phonons generated in a particle interaction 
are collected by a superconducting phase-transition ther- 
mometer [18] which is stabilized in its transition to the 
superconducting state by means of an additional heater 
structure on the thermometer [17]. By this configuration, 
the resulting change of resistance following a particle in- 
teraction can be read out with a SQUID-based readout 
scheme [19, 20] and is used for a calorimetric measure- 
ment of the interaction energy. 

An additional cryogenic detector [21] in the vicinity of 
the crystal is used to detect the scintillation light, typi- 
cally about 1% of the total energy for electrons [22]. The 
light detector response is calibrated and linearized inde- 
pendently of the scintillation light using externally cre- 
ated pulses injected into a heater structure on the ther- 
mometer [17]. Hence, each particle interaction results in 
two parameters: the energy E of the interaction as mea- 
sured directly with the calorimeter, and the light yield -g, 
which is defined as the ratio of detected light L over the 
energy E of the interaction. We measure -g in units of 
keVee/keV, where keV ee (keV electron-equivalent) is de- 
fined via a 57 Co calibration source such that E — 122 keV 
gamma rays produce photoelectrons with a light yield of 
§ = IkeVcc/keV. 

Relevant in the search for dark matter is the possibility 
to discriminate electron or gamma-induced events from 
nuclear recoils based on a strongly reduced light yield of 
the latter. For example, recoiling alphas have 18% of the 
electron light yield. Due to coherence effects, dark mat- 
ter induced events are expected to mainly cause recoiling 
tungsten nuclei, which show less than 2.5% of the elec- 
tron light yield [4, 23]. The effects discussed in this note 
are more subtle and concern slight variations within the 
population of electron and gamma-induced events. 



III. SCINTILLATOR NON-PROPORTIONALITY 

The data presented here was taken over a period of 
about one month while no calibration source was present. 
The accumulated data is shown in the light yield-energy 
plane -g (E) in figure 2 for interaction energies E up to 
200 keV. The events are mainly from a continuous elec- 
tron background due to an internal contamination of the 
crystal with 90 Sr and other beta emitters [16]. In ad- 
dition, a few gamma lines of both internal and external 
origins are superimposed; we will come back to these in 
section IV. 

Also shown in figure 2 is the mean of the light yield in 
each 2 keV energy bin, which clearly decreases with de- 
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creasing energy. This scintillator non-proportionality has 
previously been measured for CaWC>4 with small crys- 
tals under external gamma irradiation [10] at room and 
liquid nitrogen temperatures. Measurements with other 
tungstates at room temperature also show the same ef- 
fect [24, 25]. 
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FIG. 2: Electron and gamma events observed with one crys- 
tal after an exposure of 12.31 kgd (detector DAlSY/run 27). 
Number of entries per bin are according to the grey scale on 
the right. The white data points mark the mean of the light 
yield in each 2keV energy bin together with the error bars 
of the mean and the width of the energy bin; the scintillator 
non-proportionality is clearly visible. The solid (blue) line is 
a fit according to the model described in the text. 

This scintillator non-proportionality can be explained 
simply by the stopping power per unit path length 
dE/dx (E) for electrons in CaWC>4, together with Birks' 
law (equation 1). Integrating Birks' law allows a calcu- 
lation of the light yield -g (E): 

-(^)=ir^^=ir^f^v i d^(2) 

E y ' E J dE' E J dx V dx J w 

For dE/dx (E), we use the data from the estar 
database [5] shown in figure 1. With dL/dx from equa- 
tion 1, we can perform a fit of equation 2 to our data, 
which yields the parameters A and ks- This fit is also 
shown in figure 2. When the errors stated by the estar 
database are taken into account, the variation of the fit is 
within the thickness of the line in figure 2, even when al- 
lowing for a ±20% variation for energies < lOkeV. The 
model can be seen to describe the observed scintillator 
non-proportionality well. 

From the fit we obtain the Birks constants as A = 
(1.096 ± 0.003) keVee/keV and k B = (18.5 ± 0.7) nm/keV 
for electrons in our CaWC>4 crystal. Since A describes 
a linear dependence between dL/dx and dE/dx, as ap- 
propriate for higher energies, the value A ss lkeV 00 /keV 
was to be expected by construction of the keV ee unit. 
The parameter kg characterizes at what dE/dx the light 
yield g begins to saturate, and l/feg may be thought 



of as a quantity (d£/da;) saturation . Dividing by the 
density of the CaWC>4 crystal of 6.134 g/cm 3 , we have 
(d£/dx) saturation = (88 ± 3) MeVcm 2 /g. This is consis- 
tent with the range 8MeVcm 2 /g < (d-E/da;) satm . ation < 
160MeVcm 2 /g which was estimated in [4] based on a 
different method. 



IV. GAMMA QUENCHING 

The absorption of a gamma is typically a complicated 
process, where an electron is ejected from its shell, leav- 
ing a vacancy. The electron gives rise to scintillation, but 
has less energy than the incoming gamma. The binding 
energy of the electron will be released in the form of an 
Auger electron (mainly for light elements) or X-rays (for 
heavy elements), which can in turn eject other electrons, 
leading to an electron cascade for each gamma traversing 
the scintillator. In this way, the final energy transfer to 
the material following a gamma interaction takes place 
through a number of low-energy electrons, further en- 
hancing the role of the last few keV of electron tracks [7] 
(compare figure 1). 

Hence, in the presence of scintillator non- 
proportionality, we can expect a reduced light yield 
for low-energy gamma events with respect to electron 
events of the same energy, since the initial energy of 
the gamma is distributed over many electrons. Given a 
parametrization of the scintillator non-proportionality 
for electrons, the light yield for gammas has been 
calculated in the framework of this model for a few 
materials [7], although a calculation for CaWCU is not 
known to us. It thus appears possible that a non-linear 
gamma response of the light yield can be understood via 
the saturation effect of equation 1 . 

Here, we demonstrate the presence of a reduced light 
yield for gamma events in CaW04 crystals. It is already 
visible in figure 2 that the 46.5 keV line has a some- 
what reduced light yield with respect to the continuous 
electron background. In addition, at 122 keV where the 
keV C0 unit is normalized to unity with gamma events 
from a 57 Co calibration, the mean of the light yield of 
the observed electron background is (1.03 ± 0.01) keV cc , 
pointing to a reduced light yield of the gamma events in 
the calibration. 

A quantitative analysis needs to disentangle the light 
yield of gamma events from that of the electron back- 
ground. To this end, figure 3 shows the spectrum of the 
data shown in figure 2, in which three gamma lines can 
be identified. One at 46.5 keV is due to an external con- 
tamination with 210 Pb. Two more lines at 65.4 keV and 
73.7 keV are due to tungsten activated by cosmic radia- 
tion. The remaining spectrum is rather flat and domi- 
nated by the beta spectrum of 90 Sr [16]. 

We proceed as follows: for each of the observed gamma 
lines at energy E 7 , a Gaussian is fitted to the light yield 
distribution of off-peak events, which are taken from the 
intervals [E 7 - 5 keV, E 1 - 1 keV] and [E y + 1 keV, E 1 + 
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FIG. 3: Spectrum of events observed in figure 2. The flat 
background is due to the beta spectrum of 90 Sr. The three 
observed lines are labeled according to their origin: one is due 
to an external contamination with 210 Pb, the other two are 
due to cosmogenic activation of the tungsten. 



5keV]. A second population is taken on-peak, from 
[Ey— 0.4keV, £' 7 +0.4keV], which corresponds to a 2a in- 
terval around the line, given the energy resolution of the 
calorimeter at these energies. The sum of two Gaussians 
is fitted to the on-peak distribution: one Gaussian is com- 
pletely fixed to have the same mean light yield y and the 
same width a as the off-peak Gaussian, and its amplitude 
A is scaled to the smaller on-peak energy bin width. For 
the second Gaussian all parameters are left free. Fig- 
ure 4 shows the off- and on-peak distributions and the 
fits for the 210 Pb line at 46.5 keV. Clearly, the addi- 
tional gamma events from the 210 Pb line show a reduced 
light yield compared to electrons of the same energy of 
only (0.89±0.02) keV ce /keV/(0.991±0.005) keV cc /keV = 
(0.90 ± 0.02) of the electron yield. 

Since 46.5 keV gammas are absorbed within the first 
millimeter in CaW04, an alternative explanation of the 
reduced gamma light yield could be based on the follow- 
ing argument: the crystals are cut in cylindrical shape to 
have the largest crystals possible in the search for dark 
matter. Hence, events happening close to the mantle sur- 
face could show a reduced light yield due to geometrical 
light trapping in the crystal, following total internal re- 
flections on the mantle surface, in which case the 210 Pb 
line could show a reduced light yield simply due to light 
trapping. However, this is not the case, as we demon- 
strate considering sources of photons which are uniformly 
distributed throughout the crystal. 

Exposition to cosmic radiation of the crystals or the 
raw materials from which they are grown leads to an acti- 
vation of the tungsten. We observe products from the re- 
actions 182 W(p, a) 179 Ta and 183 W(p, t) 181 W. The result- 
ing isotopes will be distributed homogeneously through- 
out the whole crystal, independently of whether the ac- 
tivation takes place before or after crystal growth. 

The isotope 179 Ta decays with a half- life of 1.8 years 
in an electron capture process into 179 Hf [26]. In most 
cases, this results in a vacancy in the K-shell, which 
is quickly filled, eventually leading to a cascade of low- 
energy X-rays. In our calorimetric measurement we ob- 
serve the full binding energy of a hafnium K-shell elec- 




0.2 0.4 0.6 0. 



1.4 1.6 1.8 
Lielit Yield / keV /keV 



FIG. 4: Upper histogram: fit of a Gaussian to the off-peak 
events (only electron events) around the 210 Pb peak. Re- 
sulting fit parameters are y g = (0.991 ± 0.005) keV ee /keV 
for the light yield (in agreement with our model which gives 
0.994 keV cc /keV), and a oS = (0.116 ± 0.004) keV ee /keV. 
Lower histogram: fit of the sum of two Gaussians (full line) 
to the 210 Pb on-peak distribution. One Gaussian (dashed) 
has all parameters taken from the off-peak distribution and 
represents the electron events. All parameters of the second 
Gaussian (dotted) which accounts for the gamma events are 
left free, and the fit finds y on = (0.89 ± 0.02) keV oc /keV and 
a on = (0.09 ± 0.01) keV cc /keV. 



tron, namely 65.4 keV. The corresponding off- and on- 
peak distributions together with the Gaussian fits are 
shown in figure 5. The resulting gamma light yield 
is (0.92 ± 0.01) keV oe /keV/(1.009 ± 0.004) keV 0C /keV = 
(0.91 ±0.01) of the electron light yield, just as in the case 
of 210 Pb events. Thus, light trapping can be excluded as 
explanation for the reduced gamma light yield. 

A consistent observation can be made with a second 
line from cosmic activation: the isotope 181 W decays 
with a half- life of 121 days in an electron capture pro- 
cess into 181 Ta under emission of an additional 6.2 keV 
gamma [26]. We observe the line at an energy corre- 
sponding to the binding energy of a tantalum K-shell 
electron (67.4 keV) plus the energy of the gamma, namely 
73.7 keV. The light yield of this line is derived as above, 
and again the light yield of gammas is found to be re- 
duced to 0.88 ± 0.01 of the light yield for electron events 
of the same energy. 



V. CONCLUSION 

We have presented measurements of scintillator non- 
proportionality for low-energy electrons in CaW04. The 
crystal under investigation is cooled to mK temperatures, 
and the non-thermal phonon signal is measured with a su- 
perconducting phase-transition thermometer for a calori- 
metric measurement of the interaction energy. The light 
yield is measured with a separate cryogenic light detec- 
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FIG. 5: The Gaussian fit to the off-peak events around the 
179 Ta peak (upper histogram, oniy electron events) gives 
y oS = (1.009 ± 0.004) keV ee /keV for the light yield (in agree- 
ment with our model which gives 1.011 keV cc /keV) and <t b = 
(0.093 ± 0.003) keVe C /keV. Fitting the sum of two Gaus- 
sians to the on-peak distribution (lower histogram), one Gaus- 
sian representing electron events (dashed), the other gamma 
events (dotted), gives y on = (0.92 ± 0.01) keV cc /keV, and 
cr on = (0.034 ± 0.008) keVee/keV for the light yield of gamma 
events. 



tor. The observed scintillator non-proportionality for 
low-energy electrons is well described by a simple model 
based on Birks' law and the stopping power dE/dx of 
electrons in CaWC>4. Gamma events show a light yield 
that is reduced to ~ 90% of the light yield of electron 
events with the same energy. This is understood within 
the same model, where gamma events in a scintillator 
cause a cascade of low-energy electrons, which in turn 
produce the scintillation light. Our findings help to iden- 
tify nuclear recoil signals in experiments that use CaWO^ 
to search for dark matter. 
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